PROCEEDINGS 

°F SCIENCE 



o 

> 
O 

:z; 

(N 
(N 



00 

o 



Bound H-dibaryon from Full QCD Simulations 
on the Lattice 



Takashi Inoue* 

Nihon University, College ofBioresource Sciences, Fujisawa 252-0880, Japan 



E-mail: inoue . takashi@nihon-u . ac . jp 



for HAL QCD Collaboration 



Hadrons to Atomic nuclei 

H 




from Lattice QCD 



Using a new method recently proposed by the HAL QCD collaboration to study hadron interac- 
tions in lattice QCD, we investigate a possibility for an existence of the //-dibaryon in the flavor 
SU{3) symmetric world. A potential for the flavor-singlet baryon-baryon channel is derived from 
the Nambu-Bethe-Salpeter wave function, and the elusive //-dibaryon is shown to be a bound 
state, with the binding energy of 20-50 MeV for the pseudo-scalar meson masses in the range 
fliat 469-1 171 Me V. 



The XXIX International Symposium on Lattice Field Theory - Lattice 2011 
July 10-16, 2011 

Squaw Valley, Lake Tahoe, California 



* Speaker. 



© Copyright owned by the author(s) under the terms of the Creative Commons Attribution-NonCommercial-ShareAlike Licence. 



http://pos.sissa.it/ 



Bound H-dibaryon from Full QCD Simulations on the Lattice 



Takashi Inoue 



1. Introduction 

Lattice QCD now became a powerful method to study not only hadron spectra and QCD 
phase structure, but also hadronic interactions and multi-hadron systems. For example, hadron 
scatterings such as nn and A^A'^ have been studied by using the Liischer's finite volume method 
||]. Recently, the binding energy of helium nucleus was also measured directly in a quenched 
lattice QCD [Q]. In 2007, an approach to study hadron interactions on the basis of the Nambu- 
Bethe-Salpeter amplitude on the lattice was proposed|^, ^ : This approach (the HAL QCD method) 
has been successfully applied to baryon-baryon interactions (such as NN and YN) and the meson- 
baryon interaction (such as ^A'^)[Q, ^. 

In any of the above approaches, the spatial lattice size L should be large enough to accom- 
modate multi-hadrons inside the lattice volume. Once L becomes large, however, energy levels of 
hadrons corresponding to the scattering state in the lattice box become dense, so that isolation of the 
ground state from the excited states becomes very difficult unless unrealistically large imaginary- 
time t is employed. Recently, we proposed a new technique to resolve this issue by generalizing 
the original HAL QCD method 0, |ll]]. With our new method, information of hadron interactions 



can be obtained even without separating energy eigenstates on the lattice. 

An interesting application of the above new method is the elusive //-dibaryon which has 
been know to receive large finite volume effect on the lattice [|l^]. The //-dibaryon, predicted 



by R. L. Jaffe in 34 years ago [|13|], is one of most famous candidates of exotic-hadron. The pre- 
diction was based on the observations that the Pauli exclusion principle can be completely avoided 
due to the ^2L\ox-sing\&t{uuddss) nature of //-dibaryon, together with the large attraction from one- 



gluon-exchange interaction between quarks 1 13, 14]. Search for the //-dibaryon is one of the most 



challenging theoretical and experimental problems in the physics of strong interaction. Although 
deeply bound //-dibaryon with the binding energy Bh > 7 MeV from the AA threshold has been 



ruled out by the discovery of the double A hypernuclei, ^He [15|, there still remains a possibility 



of a shallow bound state or a resonance in this channel 1 16]. 

In this report, we first explain our new method to resolve the difficulty of the finite volume 
effect. By using the method, we then study the //-dibaryon in lattice QCD. To avoid unnecessary 
complications, we consider the flavor SU (3) hmit in this study. We find a bound //-dibaryon with 
the binding energy of 20-50 MeV for the pseudo-scalar meson mass of 469-1 171 MeV. Compar- 



ison of our result with those by other groups [ ]17| , |1 8| ] as well as possible //-dibaryon with flavor 
SU (3) breaking are also discussed. 



2. Formalism 

In the original works [Q, ^, the Nambu-Bethe-Salpeter (NBS) wave function for the two- 
baryon system with energy E, 

Mr,t)='£{0\Bi{x + r,t)Bj{x,t)\B = 2,E), (2.1) 

X 

has been employed to study baryon-baryon (BB) interactions. A correlation function *F(r,f) for 
two baryons can be expressed in terms of ^E{r,t) in a finite volume as 

^{r,t) = Ag,^E,Xr)e-'''^'' + Ais,fe,(r)^'-^'»'^ • • • (2.2) 



2 



Bound H-dibaryon from Full QCD Simulations on the Lattice 



Takashi Inoue 




Figure 1: Example of Nambu-Bethe-Salpeter wave function of two-baryon system measured on the lattice 
with L ~ 4 fm. Left (Right) panel shows the function for two-baryon in the flavor 27-plet (singlet) channel. 
All data are normalized to unity at the origin. 



where Eg^ and Sist are energy of the ground state and the first excited state, respectively, andAgr and 
Aist are the corresponding coefficients. Hereafter we call *F(r,f) a time(f)-dependent NBS wave 
function. In principle, ^ir,t) is saturated by the ground state contribution for {E\^i — £'gr) x f » 1, 
so that the wave function of the ground state can be extracted. If the ground state is a scattering 
state or a weakly bound state, the energy differences Ei — Ej are relatively large (several hundred 
MeV) in a small volume (e.g. L ~ 2 fm), and hence the condition can be fulfilled at relatively small 
? ( f ~ 1.2 fm) which we can access in actual lattice simulations. 

In a large volume, however, much larger t is required. Fig. [I] shows the f-dependent NBS wave 
function of two-baryon measured in our lattice QCD simulation on L ~ 4 fm lattice, at sink-time t 
around 10 in lattice unit (a ~ 0. 12 fm). Shapes of the function change as the sink- time t increases 
from 10 to 12, showing that is not saturated by the ground state contribution in this range 

of t. Since the energy difference tends to decrease as L^^( with assumption that a deeply bound 
state is absent), t around 40 in lattice unit (f ~ 5 fm in physical unit) may be necessary for the 
ground state saturation. We need huge statistics to extract signals at such large t, however, since 
the signal to noise ration of *P(r,?), which includes 4 baryon operators, becomes bad as t increases. 
Such a calculation would be unacceptably expensive even with today's computational resources. 
We therefore conclude that it is practically impossible to achieve the ground state saturation for the 
two-baryon system in large volume, unless we employ some sophisticated techniques such as the 
variational method with optimized sources. 

To overcome this difficulty in large volume, the HAL QCD collaboration has recently proposed 
alternative method [|^, 11 1 as summarized below. Within the non-relativistic approximation, wave 
functions we consider satisfy the Schrodinger equation with non-local but energy-independent po- 
tential U {f,r'). For the lowest two energy eigenstates, it reads 



V2' 

M1+M2- — 

2/x 

M1+M2- — 



£ist0ist(r)e 



(2.3) 
(2.4) 



where Mi 2 represent masses of two baryons and /i is the reduced mass of the two baryons. Since 
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Figure 2: Example of baryon-baryon potential extracted from QCD in our new method. Left (Right) panel 
shows the potential for two-baryon in the flavor 27-plet (singlet) channel. One can see independence of the 
potential on the sink-time t and the volume in the left and the right panel, respectively. 



these equations (and those for other E) are linear in <^e, ^ij-,t) = Ln^«fe(^)^ satisfies 

d 



M1+M2- — 



^{r,t) + jcP^U{r,f)^>{r',t) = -^"^{^^t). (2.5) 



Using this equation, we can extract U{r,f') from the f-dependent NBS wave function *F(r,f) at 
the moderate value of t, which can be easily calculated in lattice QCD simulations. This is our 
new method to study hadron interactions in lattice QCD without isolating energy eigenstates. In 
practice, we expand the non-local potential U {r,r') in terms of velocity such that U {r,r') = {V{r) + 
0(V))5^(r, — r'). In this paper, we only consider the leading term V{r) of the velocity expansion. 

Fig. § shows examples of baryon-baryon potentials extracted by the new method. In the 
left panel, we show the spin-singlet and flavor 27-plet potential V^^'^\r), which is found to be 
independent of the sink-time t. (Note that V{r) becomes a function of r = [r| for the spin-singlet 



sector.) This result shows not only a success of the new method using eq.([2.5|) but also a goodness 
of the leading order approximation of the velocity expansion, since the t independent potential 
V{r) can be extracted from the manifestly t dependent function In the right panel, the 

flavor singlet potential V^^^ (r) is plotted on L = 2, 3 and 4 [fm] lattices. Since the potential remains 
almost identical between L = 3 and 4 fm, the potential obtained on L = 4 fm is considered to be 
volume independent. This observation is consistent with a fact that the interaction range of the 
potential, which is around 1.5 fm from the figure, is smaller than the half of the lattice extension 
L = 4 fm. Once we obtain a volume independent potential, any observable of the system, such as 
binding energy and scattering phase shifts, can be extracted by solving the Schrodinger equation in 
the infinite volume. 



3. Lattice QCD setup 

For dynamical lattice QCD simulations in the flavor SU (3) limit, we have generated ensembles 
of gauge configurations on a 32^ x 32 lattice with the renormalization group improved Iwasaki 
gauge action at j3 = 1.83 and the non-perturbatively 0{a) improved Wilson quark action at five 
different value of quark mass. The lattice spacing a is found to be 0.121(2) fm and hence lattice 
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Table 1: Summary of lattice QCD parameters and hadron masses. See the main text for details. 



a [fm] 


L[fm] 


^uds 


Mp.s. [MeV] 


Mvec [MeV] 


Mbar [MeV] 


A^cfg/A^traj 






0.13660 


1170.9(7) 


1510.4(0.9) 


2274(2) 


420/4200 






0.13710 


1015.2(6) 


1360.6(1.1) 


2031(2) 


360/3600 


0.121(2) 


3.87 


0.13760 


836.5(5) 


1188.9(0.9) 


1749(1) 


480/4800 






0.13800 


672.3(6) 


1027.6(1.0) 


1484(2) 


360/3600 






0.13840 


468.6(7) 


829.2(1.5) 


1161(2) 


720/3600 



size L is 3.87 fm. Hadron masses on each ensemble are given in Table together with other 
parameters such as the quark hopping parameter Kuds, number of thermalized trajectory A'traj and 
number of configuration A^cfg- 

On each gauge configurations, the baryon two-point and four-point correlation functions are 
constructed from quark propagators for the wall source with the Dirichlet boundary condition in 
the temporal direction. Baryon operators at source are combined to generate the two-baryon state 
in a definite flavor irreducible representation, while the local octet-baryon operators are used at 
sink. To enhance signal, 16 measurements are made for each configuration, together with the 
average between forward and backward propagation in time. Statistical errors are estimated by the 
Jackknife method with bin size equal to 12 for the K^js = 0.13840 and 6 for others. 

4. Results 

We here consider the system with baryon-number B = 2 in the flavor-singlet and = 0+ 
channel, to investigate //-dibaryon. The left panel in Fig. ^ shows the baryon-baryon potential 
in this channel at five values of quark mass: it is entirely attractive and has an "attractive core" 
[^]. This result is consistent with the prediction by Jaffe and also by the quark-model in this 
channel. The figure also indicates that the attractive interaction becomes stronger as the quark 
mass decreases. 

By solving the Schrodinger equation with this potential, we have found one bound state in this 
channel ^n\\. The right panel in Fig. ^ gives energy and size of this bound state, showing that a 
stable //-dibaryon exists at this range of the quark mass in the flavor SU (3) symmetric world. 

As shown in the right panel of Fig. [l], the f-dependent NBS wave function for this channel goes 
to non-zero value at large distance, contrary to a naive expectation for bound state wave functions. 
The non-vanishing wave function at large distance is due to the contribution from excited states; 
they correspond to scattering states in the infinite volume and do not vanish at large distance. In 
other words, the ^-dependent NBS wave function is a superposition of the bound state and scattering 
states. Even in such a case, our new method is shown to work well. 

The binding energy Bh = —Eq for the //-dibaryon ranges from 20 to 50 MeV and decreases 
as the quark mass decreases. The rms distance yir^ is a measure of a "size" of the //-dibaryon, 
which may be compared to the rms distance of the deuteron in nature, 1.9 x 2 = 3.8 fm. Although 
quark masses are different between the two, this comparison suggests that //-dibaryon is much 
more compact than the deuteron. 
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Figure 3: Left: Baryon-baryon potential in the flavor singlet channel. Right: The ground state of the system. 
The energy £0 is measured from two-baryon threshold. The bars indicate statistical error only. 
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By including a small systematic error caused by the choice of sink-time t in the ^-dependent 
NBS wave function, the final result for the //-dibaryon binding energy become 

(4.1) 
(4.2) 
(4.3) 
(4.4) 
(4.5) 

with the statistical error (first) and the systematic error (second). A bound //-dibaryon is also 
reported by the full QCD simulation with a different approach from ours [|l^]. The obtained bind- 
ing energy from the AA threshold Bh = 16.6(2. 1)(4.6) MeV at {Mjc,Mk) ^ (389,544) MeV is 
consistent with the present result. 

A deeply bound //-dibaryon is ruled out by the discovery of the double A hypernucleus. The 
binding energy Bh in this paper, however, should be interpreted as the binding energy from the BB 
threshold averaged in the {S,I) = (—2,0) sector. In the real world, the AA, NE and ZZ thresholds 
in this sector largely split due to the flavor SU (3) breaking. We therefore expect that the binding 
energy of the //-dibaryon measured from the AA threshold in nature becomes much smaller than 
the present value or even the //-dibaryon goes above the AA threshold. Our trial calculation using a 
phenomenological flavor SU (3) breaking indicates that the //-dibaryon state becomes a resonance 
state above the AA threshold but below the NE threshold, as the flavor SU (3) breaking reaches 
the physical value. To make a definite conclusion on this point, however, we need lattice QCD 
simulations at the physical point, together with the coupled channel analysis. We have already 
developed the formula for this purpose |19] and tested the method in numerical simulations [pO|]. 
A study towards this final goal is now in progress. 
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